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1.0 Introduction 


Since 1980, personnel of Applied Research, Inc. have suppor- 
ted NASA at Marshall Space Flight Center in the measurement of 
the atmospheric backscatter coefficient, 6, with an airborne CO 
Laser Doppler Velocimeter (LDV) system operating in a continuous 
wave, focussed mode. A method, called the Single Particle Mode 
( SPM ) algorithm, has been developed from concept through analysis 
of an extensive amount of data obtained with the system aboard a 
NASA aircraft. The SPM algorithm is intended to be employed in 
situations where one particle at a time appears in the sensitive 
volume of the LDV. In addition to giving the backscatter coeffi- 
cient, the SPM algorithm also produces as intermediate results 
the aerosol density and the aerosol backscatter cross-section 
distribution. 

A second method, which measures only the atmospheric back- 
scatter coefficient, is called the Volume Mode (VM) and was 
simultaneously employed in obtaining the aforementioned data. The 
results of these two methods generally differed by slightly less 
than an order of magnitude. The purpose of this report is to 
examine the measurement uncertainties or other errors in the 
results of the two methods. 

A review of the basis of each method is given in Section 2, 
including a discussion of the computer programming implementation 
of the SPM. A discussion of error inherent in each is given in 
Section 3, with conclusions summarized in Section 4. 
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2.0 Review of the Algorithms for Determining the Backscatter 

Coefficient 

For convenience, the basis of the VM and SPM algorithms for 
obtaining the atmospheric backscatter coefficient B is presented 
here. It will be seen that the VM method works under more 
general conditions than the SPM method, but that the latter gives 
more information, namely the aerosol density and backscatter dis- 
tribution. Each method has its own special calibration 
requirements . 

2.1 Volume Mode Algorithm 

The VM algorithm will now be derived from the response of 
the LDV to a single particle in its sensitive volume. This 
expression will be compared to the well known result for the 
signal-to-noise, S, of a focussed, cw LDV operating in the con- 
ventional "volume mode" with very many particles in its sensitive 
volume: S = B G v where 

3 = y* a n(a)d a 

a = single particle backscatter cross-section (m 2 ) 
n(o) = backscatter cross-section distribution (m -5 ) 

Gy = volume mode gain factor determined by 
calibration (m). 

In this derivation, the system will be assumed to be aircraft 
borne with optical axis perpendicular to the aircraft velocity 
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vector. In this case each aerosol particle passes through the 
sensitive volume perpendicular to the optic axis with velocity v. 
See Figure 1. A particle at x,y,z with cross-section o produces 
a signal-to-noise at time t given by 

S( 0,t) = ou(x,y,z) 

where u(x,y,z) defines the LDV sensitive volume. That is, 
requiring S>1 defines a sensitive volume dependent on o . Extend- 
ing this expression to many particles with cross-section a moving 
in the positive x direction, each with initial position x Qi such 
that 

x i = vt + x oi gives 
S(a,t) - oE u ( vt + x oi ,y i ,z i ) 

Further extending to a time average over M time intervals At 
gives 

M 

S( 0 ) « ^S( o,t . )/M 

j=l 

= a ^u(y i ,z i )m i /M 
i 

with u(y i ,z i ) = ^ u(vt^ +x Qi ,y i ,z i )/m i 

j 

and mi is the number of integration points within the sensitive 
volume. This quantity depends on the trajectory of the particle 
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y is up 

z is laser beam 

x is along aircraft velocity vector 

Particles penetrate perpendicular to the yz plane 


Figure 1. Geometry and Coordinates 



defined by the constants and z^ through the sensitive volume 
and is independent of velocity. Actually 

mi = Ax(yi ,Zi ) /v At 


where x(yi,Zi) is the width of the sensitive volume at yi and z ^ • 
Once again extending the expression by replacing the summation 
over particle number by integration over volume and the uniform 
density n( o) gives 


S(a) = n(a) ay' u (y,z)Ax(y,z)dydz 


= n(a ) a G(a ) . 


A possible dependence of the gain factor G on a is emphasized 
because of its importance in case this method is used when a 
single particle at a time is sensed. That is, if an effective 
threshold exists at S=l, the sensitive volume is defined by 
S>1, and the single particle sensitive volume is smaller than the 
many particle sensitive volume because many particles outside the 
single particle volume can build up a signal S>1. In order to 
include all particle sizes, the expression must be integrated 
over o : 

S = J' n(a ) a G(a )d a . 

Putting this in the form given at the beginning of this section 
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such that 6 can be isolated requires defining 


so that 


VM 




/• 


) o G ( o ) d o/ I n ( a ) ado 


S - 3 G 


Note that G^j<G v , the conventional many particle volume 
mode gain, if G(a) is dependent on a . This means that, 
if 6 is calculated from a measured signal S obtained from an 
average of single particle signals with 3 = S/Gy , too small a 
result would be obtained. 

The basis of the VM method for measuring the backscatter 
cross-section has been developed. However it has been shown that 
the gain factor may depend on whether data is taken in a many 
particle or single particle situation. 
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2.2 Single Particle Mode Algorithm 


The SPM algorithm functions by recording the peak signal 
from each particle which transits the sensitive volume. The 
geometry is the same as described previously. Single particle 
transits are assumed not to overlap in time. From the statistics 
of the peak signal distribution and knowledge of sensitivity con- 
tours within the sensitive volume, one may derive the particle 
density, backscatter cross-section distribution,' and atmospheric 
backscatter coefficient. The algorithm will now be described 
using discrete mathematics with direct application to the compu- 
tational programs. 

Consider an interval of the single particle backscatter 

cross-section axis from a T to a„ which covers all particles seen, 

L H 

and which itself is divided into M intervals. The backscatter 

coefficient will be taken as 

M 

B = 2> j0j 
j = l 

where 


nj = number of particles per unit 

volume within the jth interval 

Oj = cross-section at the center 
of the jth interval. 

Suppose that the total number of particles seen per unit volume 
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is D; then 


w j = n j/^ 


is a probability distribution since 


M M 

= X>-j/D = 1 

j=i : =1 


Since each Oj has a particular amount of the sensitive volume Vj 
in which the particle would give a signal above threshold S L , 
the total number of particles seen, N, is 


M M 

n = y^ n j v j = d yy-jVj 

j=i j=i 


and the density D is 

M 

D = N / ^ w j V j 

j = l 


which says that the relevant volume is an average volume weighted 
by the cross-sectional distribution Wj . Also 


M M 

3 = ^ njCj/4IT = D^ WjOj /4IT 
j = l D = 1 

= d <o>/4n 
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so that 3 is given by the particle density times the average 
cross-section. The SPM algorithm determines the probability 
distribution Wj through the relationship between the peak signal 
distribution and the sensitivity contours within the sensitive 
volume. This same information also determines the density. 

In order to understand the SPM algorithm, the concept of the 
volume Vj sensitive to a particle of cross-section Oj must be 
clear. For the geometry previously described with the laser axis 
perpendicular to the aircraft vector, this volume is a cylinder 
with axis in the direction of the aircraft velocity vector and 
length vt, where v is the aircraft speed and t is the time of 
observation. This cylinder intersects a vertical plane through 
the optic axis in an area Aj . The definition of the volume Vj 
is that particles Vj passing through the area Aj will give a peak 
signal S (signal-to-noise ) above threshold S L . The gain factor 
g(y,z) in the equation S L <^g(y,z) determines the area Aj , 
where y,z are coordinates in the vertical plane of the sensitive 
volume. Furthermore, the notation g(A) is used to mean the gain 
on the contour surrounding the area A. These concepts and the 
relationship between the a and S distributions are shown in 
Figure 2a. The line g(0) represents the maximum gain of the 
sensitive volume because the maximum gain occurs at a point. 

Three a and S intervals are shown in this figure, with boundaries 
related by S=og(0). Any other gain is represented by another 
line with slope g(a) < g(0). Notice that a particle in interval Ao 2 
may contribute to AS X and AS ^ > depending on the gain, but not 
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Figure 2a. 



Figure 2b. 



Figure 2. Peak Signal versus Backscatter 
Cross-section Diagrams 
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to AS 3 since this would require a gain larger than g(0). 

In order to relate the a and S distributions, consider 
Figure 2b. Particles with a number per unit area density n(c?) 
at o contribute signals AN within 6 S according to AN = 5An(o) 

To determine the contribution of a macroscopic interval A a to a 
macroscopic interval AS, as shown in Figure 2c, the intergral 

AN = J' r\ ( a) AA ( a) 6a 
A a 

is evaluated where n(a) is the number of particles per unit vol- 
ume per a increment at a , AA(o) = A 2 - A* , and 6o is a sub- 
interval in a . Notice that the areas can be considered 
functions of the gain, and therefore A 2 = A(S 2 /o) and 
A x = AfSj/a). This function has been determined by calibration 
and is included in the computational program (Section 2.3, 2.4 
and Appendix in tabular form in subroutine DATANAL. The method of 

obtaining these data is described in Reference 1. In order to 
perform the above integral, the probabilities behind the n(o) are 
assumed (w^ in previous discussion) and the probability of 
obtaining AN particles is calculated in subroutine OPTIM. The 
actual numbers of particles are obtained by multiplying by the 
total number of particles observed. Contributions to each AS 
interval from each possible Aa interval are evaluated. The 
assumed distribution n(a) is then varied until the best fit to 
the peak signal distribution is obtained in the least squares 
sense. With this "best fit" distribution the density, average 
cross-section <o>, and backscatter coefficient 8 are determined. 


11 



2.3 Computer Implementation 

Implementation of the 6 prediction algorithms involved 
setting several parameters and making certain assumptions. These 
parameters/assumptions have little precedence, and therefore 
Applied Research will identify these variables/assumptions as 
clearly as possible. 

1. The number of signal bins = 6. Six bins were processed 
and used for the inversion. A seventh bin was 
identified and used to store all large particle signals, 
but not considered in the inversion. 

2. The number of processor bins which made up each of the 
six signal bins were: 

Signal bin 1 = processor bins 5-8 

Signal bin 2 = processor bins 9-16 

Signal bin 3 = processor bins 17-32 

Signal bin 4 = processor bins 33- 64 

Signal bin 5 = processor bins 65-128 

Signal bin 6 = processor bins 129-255 

Signal bin 7 = processor bin 256 

Note that signal bins 1-6 contain progressively twice as 
many processor bins as the previous signal bin. This 
type of consolidation was done because the lower order 
processor bins get many more signals/bin than the higher 
order bins. 

Note also that processor bins 1-4 are not used because 
the threshold was set at 4. 
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3. The signal processor threshold is equal to 4. 

4. The volume channel bandwidth is 860 Khz. 

5. The single particle processor bandwidth is 1.5 Mhz. 

6. The LDV is always focussed at 10 meters. 

7. The number of sigma bins = 6. 

The above mentioned parameters/assumptions were employed 
throughout the data prediction results. 

In addition to the pre-set parameters, several other pieces 
of information are very important to the data predictions. This 
information is all related to calibration of either the single 
particle mode or the volume mode. The volume mode calibration 
and method of prediction will be discussed first, followed by the 
single particle calibration and method of prediction. 

For the VM calculations, two sets of calibration 
data were used. The first is identified as the "count" calibra- 
tion and this data is read in at the beginning of the main 
program and stored in the array CAL (J,I). The first index J is 
over IF gain values from 45 to 70. The second index I is over 
DBSM values from -71 to -20. The proper calibration curve was 
picked depending upon the IF gain parameter. The volume 
"count" calibration data also had to be divided by 500 in order 
to scale the counts to one average count value.* The volume 
"count" calibration data relates "counts" in the volume channel 


♦Volume mode calibration data was performed by W. Jones of NASA. 
This information was relayed verbally to Applied Research, Inc. 
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to DBSM values. The second set of calibration data required is a 
measurement of the signal response of the LDV as a function of 
longitudinal distance while focussed at 10 meters. This data is 
shown in Figure 3. This data was used to calculate I^ff at 10 
meters. (Reference 1) Longitudinal distance means along the 
optical axis of the LDV. 

The volume 6 was calculated using the following equation: 

3 = (Vs/V n )f r /(S^I feff ) 

V g - integrated volume signal 

V n - integrated volume noise 

- signal-to-noise from a sandpaper disk 
at 10 meters focus with noise scaled from 
the calibration bandwidth of 100kHz to 
860kHz, the processor bandwidth, giving 
1.84x10 6 . 

f r - bidirectional reflectance of a sandpaper 
disk equal to .016. 

L e ff ” effective length of the LDV while focussed 
at 10 meters, found to be 64.1 cm. 

L 0 fff and f r were all determined prior to beginning the data 

predictions and remain fixed throughout. The volume signal and 
noise were calculated for each data point. The LDV data had 
alternating samples of noise data throughout each of the flights. 
This was accomplished by dithering a mirror in the optical train 
thereby effectively "losing" the signal. The raw output of the 
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volume 6 channel on the LDV processor is called "counts". This 
count has to be divided by the number of times the channel was 
scanned in order to provide the average count for the time 
interval. The average count is then compared to the count 
calibration data (using the correct IF gain) to obtain a DBSM 
signal. The assumption was made that the volume channel count 
(data record) actually was signal plus noise while the noise 
particle count was noise only. Therefore to obtain the volume 
signal (Vs) the noise signal (V n ) was subtracted from the signal 
plus noise. This provided the V s and V n required to complete the 
volume 3 prediction. 

Two pieces of calibration data are required to complete 
calibration for the single particle 6 calculation. These data 
are the correlation between signal value and channel number on 
the signal processor and the gain (equals signal/noise/cross- 
section) versus area curve. 

The correlation between channel number and signal value is 
important to establishing what size particles are being seen by 
the processor. The "signals" here are tied to the calibration of 
the area curve discussed next. It is not known how this cali- 
bration data may depend upon the IF gain. All predictions were 
obtained using one calibration curve. Figure 4, which has an 
IF gain of 57. 

The gain versus area curve is crucial to making the correct 
single particle 3 predictions. This data was taken by Applied 
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Figure 4. Signal Processor Channel vs. Signal Level 


Research, Inc. and presented in Reference 1. This data resulted 
in a relative mapping of the sensitive focal volume of the LDV in 
a vertical plane containing the optical axis. "Relative" means 
that the scale of the map (signal/noise/sigma) was undefined. 

The remaining procedure then was to define the scale of the gain 
versus area map. This was done by using the following equation: 

Sd/f r = / Cg ' (A) dA 

where S d is the signal return from a sandpaper disk at 10 meters, 
is the bidirectional reflectance of a sandpaper disk, and g' is 
the relatively scaled gain data. The values used were: 

S d = 1.053 x 10 6 
f r = .016 

/g'dA = 1.44 x 10 6 

The S d value is the calibration value of 10 7 * 2 scaled from 100kHz 
to 1.5MHz. The value of C determined was : C = 45.99. This 
value is used in subroutine DATANAL in a data statement as 
variable TRANS to scale the area curve to the proper values. 

Also as a result of this determination of a scale factor C, the 
single particle gain at an area of 0 can be used from the 
equation S/a/ a _ =g(0). to relate a S/N value S to a o value; 

i.e. o= S/g(0). This relationship is used in the main program 

* 

about lines 218-220 with a value of 1.1037 x 10 7 /cm 2 for g(0). 

The result of the inversion is to return an average sensi- 
tive volume cross-sectional area (WV) and average particle 
cross-section (WS). The equation 

8 = ( ISUM/4 ID (WS/FLTL/WV ) 


18 



was used to compute the single particle B , were ISUM is the 
number of particles in the histogram. The flight length, FLTL, 
is computed by using the data recorded in each record (or 
records) for the elapsed time and flight speed. 

The results of the computer implementation were presented in 
Reference 2. The results presented were all obtained using a 
power law solution for the a distribution in the inversion 
subroutine OPTIM. This was done because early test runs 
indicated that the power law solution was nearly always the best 
fit. Since run time was becoming a problem for processing the 
many flights, the exponential and log-normal solutions were 
bypassed. This "bypass" can be eliminated by removing line 
603.07 - IF ( IP . EQ . 2 ) RETURN in the subroutine OPTIM. 

Figures 5-7 contain a flowchart of the main program and 
various subprograms which constitute the computer implementation 
of the 3 prediction algorithms. The Appendix has a short descrip- 
tion of each program of the 3 prediction algorithms and also 
contains a listing of the main program and associated 
subprograms. 
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Flew Chart of Computer Implementation 
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Figure 5b. Flew Chart of Computer 

Implementation (continued) 
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Figure 6a. Flow Chart of Computer Implementation 

of Single Particle Inversion Subprogram 
OPTIM 
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Figure 6b. Flow Chart of Computer Implementation of 
Single Particle Inversion Subprogram 
OPTIM (continued) 
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Figure 7. Flowchart of Computer Implementation of 
Single Particle Inversion Subprogram 
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2.4 Running The Data Prediction Algorithms 


The 8 prediction algorithms are stored on disk on the Sigma 
V computer in account BILBRO, EB23175. The file name is DATANAL. 

A batch file which is used with the code is stored under BXD. Two 
input files are required, they are unit numbers 5 and 12. File 
unit number 5 must be LDV processor data in the format shown in 
Figure 8. This defines the data for which predictions are 
to be made. File unit number 12 is volume 8 calibration data and 
is stored under the file name CALF1L4 in account BILBRO, EB23175. 
Units 6 and 10 are used as output unit numbers to obtain the 
printed results of the algorithm. Unit 13 is also an output unit 
but is used to store results for plotting and should therefore be 
given a file name by the user. Unit 13 produces one record for 
each 8 prediction during the flight in the following format: 

Wr i t e ( 1 3 ) LH , LM , LS , LF , S INGLEB , VOLB , WS , FLTL , ELS , W V . 

The definition of these variables is: 

LH - Hours for this prediction 
LM - Minutes for this prediction 
LS - Seconds for this prediction 
LF - Fractional seconds for this prediction 
SINGLEB - Single particle 8 prediction 
VOLB - Volume channel 8 prediction 
WS - Predicted average cross-section 
FLTL - Estimated flight length 
ELS - Number seconds for which data-recording 
occurred for this prediction 
WV - Predicted average focal "volume" cross-section 


26 



3.0 Measurement Errors 


In the application of the previously discussed VM and SPM 
algorithms to the processing of the NASA flight data, two general 
classes of error are considered. The first class, systematic 
errors, involves errors of concept, principle, or application. 

The second class concerns statistical errors or uncertainties 
which result from lack of knowledge or from fluctuating random 
variables. Systematic error is the most difficult to handle 
since even its presence is not always obvious. In the case of 
the data under consideration, an average offset of the results 
from the two algorithms of about one order of magnitude indicates 
systematic error. However, more than one such error could be 
present. To uncover these errors, a reconsideration of concepts 
and procedures is necessary. 

Statistical error is more easily treated. With estimates of 
the statistical uncertainty in independent variables, the 
uncertainty in the final result may be obtained. In the present 
case, also the number of particles processed per output 8 value 
is expected to govern statistical fluctuation in 8 . This result 
has been analyzed. 

3.1 Volume Mode Errors 

The VM algorithm has been derived in Section 2.1. Possible 
systematic errors in the backscatter coefficient calculation may 
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reside in 


a) the processor output 

b) the interpretation of the processor output 

c) the concept for the volume mode gain when 
used with single particles (as discussed in 
Section 2.1) 

d) a systematic offset in the parameters used 

to calculate the backscatter coefficient value 

Possible errors in the processor output or the curves supp- 
lied by NASA which translate output into signal level are not 
considered here, but must be investigated in the laboratory. 

The interpretation of the processor output has been 
discussed in Section 2.3 and no error has been discovered. 

As discussed in Section 2.1, if the gain factor G(o) is 
dependent on a as a result of a decreased sensitivity volume for 
smaller o , calibration of the "single particle volume mode" 
should not be treated as the "conventional volume mode." As 
mentioned in Section 2.1, G VM <G V implies that B values calculated 
with the gain factor Gy for the conventional volume mode give too 
small a value by the ratio Gvm/ g V* Evaluation of Gvm requires 
laboratory determination of the sensitive volume for each value 
of a , integration of the single particle gain over this volume, 
and then integration of this result over the backscatter cross- 
section distribution. Therefore, in measurement situations, this 
effect would require determination of the a distribution. An 
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effect very similar to this occurs in the SPM where the aerosol 
density is calculated by dividing the number of particles seen 
by the average volume, calculated by averaging the sensitive 
volume elements for the SPM over the o probability distribution. 
It is not clear whether the VM sensitive volume elements would be 
the same as those for the SPM since signal acquisition require- 
ments are different for the two modes. One has the feeling that 
this effect would contribute errors of less than an order of 
magnitude. This question might be answered with computer 
modeling. 

The backscatter coefficient is calculated from the 
expression 

3 - S/Gy 


where S is the signal to noise found from the data as explained 
in Section 2.3, and Gyis found by calibration with 


where 


G V “ ^^ff / f r 

S d = signal to noise return from 

a calibration disk of bidirec- 
tional reflectivity f r 

L eff = Jg(x,y,z)dxdydz//g(x,y,0)dxdy 
(see Figure 1) 
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Lgff = the effective length of the 
sensitive volume 

Errors or uncertainties in B may be evaluated as 

6B/3 = 6 S/S + 6S d /S d + 6L eff /L eff + 6f r /f r - 
The change AB in units of orders of magnitude is defined as 

AB = log( B + 6 B ) ~ log B 
= log ( 1 + 6B / 6) 

Each of the above error sources is considered separately. 

Error due to S: 

The signal out of the processor and corrected by noise sub- 
traction is assumed here as exact. A possible interpolation 
error of less than ldBm occurs in obtaining the dBm level from 
the volume channel counts. This yields a AB found as 

AB = log ( 1 + 6 S/S) = log(S+6S)- log S 

= f (S+6s) dBm -W /10 

= .1 
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Error due to % : 


The calibration signal error from the rotating sandpaper 
disk at the sensitive volume waist would not be likely to have a 
reading error of more than 2 dBm, from inspection of the varia- 
tion of the calibration data. Thus 


A6 = .2 


Error due to f r : 

The bidirectional reflectance of the aluminum oxide target 
for incidence and reflection at 45 degrees was taken to be 
. 016/sr. An uncertainty of ±25% has been quoted by workers in 
this field (Reference 3). This gives 

A6 = .1 

Error due to f : 

The error in the determination of I^ff is less than 10% from 
consideration of the methods used. This gives 

A6 = .04 

In summary, the possible errors due to uncertainty in the 
quantities used to calculate the backscatter signal in the volume 
mode are given in the following table: 


31 



Source 


S 


Sd 

fr 


A8 

(orders of magnitude ) 

.1 


.2 

.1 


.04 


These uncertainties may be positive or negative, and combine to 
less than half an order of magnitude at maximum. The error given 
for S considers only the uncertainty in interpolating the 
processor output. Systematic errors are not included in these 
numbers, but once identified, their effects would be calculated 
with the above formulae. 
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3.2 Single Particle Mode Errors 


As discussed in Section 2.2, the backscatter coefficient 
for the SPM may be expressed as 

3 = D<a>/ 4 II 

where D is the particle density and <o> is the average single 
particle backscatter cross-section. These quantities are deter- 
mined from the peak signal distribution and the gain contours 
versus enclosed area function resulting from mapping the 
sensitive focal volume area in the vertical plane of maximum 
gain - the yz plane of Figure 1. 

Possible systematic errors for this mode concern 

a) the processor output 

b) the interpretation of the 
processor output 

c) the single particle per transit 
time condition on the applica- 
bility of the algorithm 

d) the gain versus area function 
determined by calibration 

e) operation of the SPM algorithm 

Random errors and measurement uncertainties will be discussed 
following discussion of these topics. 

The processor output will be assumed valid, subject to 


33 



laboratory tests, except possibly for point c). The interpre- 
tation of this output is described in Section 2.3 and questions 
are raised under point d). Some implications regarding conditions 
of validity for the SPM algorithms also follow from the 
discussion of point d) below. The SPM algorithm has been found 

valid for predicting 3, subject to random errors, as previously 
documented . 

In the course of the present work, considerable effort has 
been spent in an attempt to improve the application of calibra- 
tion procedures which result in the single particle gain versus 
area function, with the result that a significant shift in the 
previously processed backscatter coefficient values is indicated. 
Originally, calibration plans for this program included a special 
calibration device which would emit single particles of known 
cross-section. Time constraints from the flight plans forced 
this effort to be interrupted and an alternate calibration 
procedure was sought involving the spinning disk. Under this 
contract this procedure has been reconsidered and improved. 

Contours of constant single particle gain, g(x,y,z), were 
mapped out in the x=0 (vertical) laser beam plane within the 
sensitive volume, in a previous contract. Only relative values 
of these contours were determined by this mapping procedure, with 
overall normalization to be determined by another method. A more 
accurate way of obtaining this normalization follows. 

The conventional volume mode gain G v is defined as the 
integral over the single particle mode gain g(x,y,z) and 
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determined by calibration with a spinning disk at the waist as 


G v - /g(x ,y , z )dxdydz 


" S d L eff /f r 


The desired integral to be evaluated for normalization of the 
gain versus area function is 


I 


= J g(o,y,z)dydz 

-/< 


g ( A ) dA 


where a change of variables transforming the integral from a two 
dimensional to a one dimensional integral is implied. The varia- 
ble A is the total area within closed contours of g. With this 
integral evaluated, g is related to its unnormalized version g' 
by 

g = ag' 


where 

a = I/Jg' ( A ) dA . 


The value I can be determined from Gy if the sensitive volume is 
assumed cylindrically symmetric about the optical axis of the 
beam. Representing the integral for G v by transforming to 
cylindrical coordinates gives 
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°v 


■ / 9 (p - 


z)pdp d0dz 


lJq(P, 


= 2n / g(P, z ) pdpdz 


= 211 p av J*g(p,z)dpdz, 


This last integral is one-half of the desired integral I. Also 


p av h (p , z )pdpdz/Jg (p , z )dpdz 


is an average radius of the sensitive volume. Hence 


^ = n Pav i = s L eff /f r 


and from 


I = a/g' (A)dA = al' 
one gets for the normalization factor 

a = tS a /(f r IM]L eff /(np av ) 

= Leff/(n P av ) 


where a Q was the normalization used in the processing of the 
flight data. The integrals for P av were numerically evaluated 
to give p av = .0265cm. (Compare this with Figure 9.21 of 
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Reference 1.) This gives 


a = [ ( 64 . lcm)/( . 0265 II cm) ]a Q 
= 770a o 

Any adjustment in the normalization of the gain versus area 
function affects the sizes of the particle distribution which can 
be detected and associated with a particular peak signal distri- 
bution. For the power law backscatter cross-section distribution 
function, which best fits the processed flight data, an adjust- 
ment "a" divides the average cross-section, but does not change 
the density. (See subroutine OPTIM). Thus 

B ■ <cr> D/4n 
= 8 o /770 

and indicates that all calculated single particle B values should 
be decreased by about 2.9 orders of magnitude. A calibration 
factor of the same order of magnitude is indicated by using the 
incoherent beam profiles of Figures 8.3 and 8.4 of Reference 1. 
This implies that the system is sensitive to cross-sections as 
small as .00166 square microns since threshold peak signal 
levels required cross-sections of 1.28 square microns at the 
original calibration. 
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The signal-to-noise for a given cross-section is unchanged 
by this rescaling since its effect on g and a cancel for a power 
law a distribution. For example, around record 6 of flight 
number 9, processing of 10,000 particles gives a signal-to-noise 
of about .14 for the threshold in this case. To obtain this 
number one must use the calibration curve supplied by NASA for IF 
gain of 57. Flight 9 had an IF gain of 63. It is clearly these 
low signal-to-noise numbers combined with the higher single 
particle gain values which are responsible for such low 0 values. 

This adjustment in the normalization of the gain means that 
the system is sensitive to smaller particles than previously 
indicated, without undergoing an adjustment in the required 
signal-to-noise threshold. The backscatter coefficient is conse- 
quently smaller because no corresponding adjustment in density 
occurs. This gives a result which is about two orders of magni- 
tude lower than the VM measurement. Nominal densities calculated 
with the SMP algorithm give about 1 particle per 100 cubic centi- 
meters, which seems rather low for distributions with an average 
cross-section of .007 square microns, as in the flight 9 case 
above. This density would be more appropriate for particles of 
greater than 1 micron diameter. Therefore, since the gain pro- 
file of the sensitive volume is fairly well established (coherent 
and incoherent methods give the same order of magnitude for the 
scaling), there is reason to question the signal-to-noise values 
assigned to the signal bins of the processor output. The method 
of calculating these has been rechecked with no mistake found. 
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The actual calibration curve provided by NASA which assigns 
signal level to a bin should be rechecked. Also the variation of 
these assignments with IF Gain should be considered. An increase 
in signal values through recalibration of the processor gain 
could bring SPM and VM results into agreement. 

The following consideration suggests another possibility. 
With the same particle distribution, the VM result predicts 
densities of about two orders of magnitude higher. For this case 
of one particle per cubic centimeter, the system may see more 
than one particle per transit time since the sensitive cross- 
sectional area is 42 square centimeters. A coincidence rejectior 
mechanism could give erroneous results since no compensation is 
made for rejected cases; or a coincidence rejection mechanism 
which is not properly responding .could interpret clusters of 
particles as single particles, thereby decreasing the density. 
Hence it may be possible that the SPM was used outside the con- 
ditions of validity without proper compensation being made. 

Random errors, or errors of uncertainty of measurement, in 
the single particle prediction can result from the gain versus 
area function determination, the interpretation of the peak 
signal bins, the operation of the algorithm, and the number of 
particles seen. As before, the uncertainty in 6 in units of 
orders of magnitude due to a component uncertainty A8 is 

AB = log (1 + 68/ 8 ) 
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As documented in an informal memo, the algorithms produce an 
uncertainty of less than 5%. 

Error due to algorithm: 

A8 = .02 

Random error due to interpolation of signal bin calibration 
is 1 dBm. (This does not include possible systematic error due 
to this calibration data discussed earlier). For small errors 
this has an approximately linear affect on the a interval, and 
therefore on <a>. Hence 

<58 = 26% 

Error due to interpolation of peak signal calibration: 


A8 = .1 

The gain versus area function was previously discussed in 
terms of systematic error. Any alteration of this function by a 
factor "a" goes directly through to <o> and then to 68 . It is 
helpful to derive an estimate for the uncertainty in "a" due to 
measurement errors by referring to the previous expression 

a = a o L eff/ p av n 

where 

a o “ ^< 3 / f r I ' 
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and I' is the integration over the unnormalized contours. The 
uncertainty in I 1 is difficult to estimate precisely but is 
within 10%. Using the previous values for and f r , and taking 
p av to have the same uncertainty as I' and L e ff gives the follow- 
ing results for random errors due to uncertainties in measure- 
ments used to calculate the backscatter coefficient in the single 
particle mode. 

A6 

Source (orders of magnitude ) 


Algorithm 

.02 

signal interpolation 

.1 

S d 

.2 

fr 

.1 

f’ef f 

.04 

Pav 

.04 

I' 

.04 


Errors due to particle number fluctuations are not treated 
in this format but considered in the following paragraph. 

Errors resulting from fluctuations due to small particle 
samples have been examined by comparing with the results from the 
case when a 10/000 particle minimum was used. For processing 
runs utilizing various sample sizes less than 10,000, a running 
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time-weighted average of the results sufficient to include 10,000 
particles was developed, and considered a "local average". This 
result was used in two ways. The difference of this local 
average from the 10,000 particle case was considered the error in 
the local average. (The local average value nearest in time to 
the 10,000 particle value was used). Also, the fluctuation from 
the local average was calculated from the differences of the 
small particle number data from the local average. This method 
was found to be reasonable if too sharp variations in the data do 
not occur. Figure 9 shows these results for a region of flight 
16 with the sharp 3 spike at time 878 removed. These results 
show rms deviations in units of orders of magnitude. They 
indicate that the local error can be maintained with .1 orders of 
magnitude with particle samples of 1000 minimum. Fluctuations of 
the "instantaneous" data from the local average do not approach 
this value until around 5000 particles minimum are processed. In 
this case the SPM has lower deviations and fluctuations than the 
VM data. (The VM data was also compared with itself in the 
manner described above). Figure 10 shows the same flight with 
the B spike included, with results which do not settle down as 
well. Figure 9 shows the deviations of the local average to be 
smaller than the fluctuations of the data from the local average. 
Figure 10 shows some of the SPM versus VM behavior. 

Figure 11 shows the same data from a region of flight 9. 

Here the fluctuations are lower than the deviation of the local 
mean, and the results do not settle down as well as the 10,000 


42 



particle case is approached. For the SPM the .1 order of magni- 
tude error is approached at 1000 particles, while 3,000 particles 
seem to be required for the VM result. 
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Figure 9. Statistical Error in Backscatter Coefficient as a 
Function of Number of Particles Processed for 
Flight 16 (Spike Removed) 
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Figure 10. Statistical Error in Backscatter Coefficient as a 
Function of Number of Particles Processed for 
Flight 16 (with Spike) 
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Figure 11. Statistical Error in Backscatter Coefficient as a 
Function of Number of Particles Processed for 
Flight 9 
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3.3 Anomalous Cases From Processed Data 


Two pathological data prediction cases can be identified 
in Reference 2 and are shown in Figures 12 and 13. These cases 
are Flights 4 and 13. Flight 4 is unusual in that initially the 
single 0 predictions were much lower than the volume 0 and then 
about halfway through the flight the situation reversed. It was 
found that this was the case because the airspeed value was 
incorrect. Since airspeed is used to calculate flight length and 
then single particle 0 , the predicted results were much too low. 
Later on in the flight the problem was corrected and the 
predictions reversed the earlier trend. 

Flight 13 also showed the trend of single 0 predictions less 
than the volume 0 predictions. A closer look at the computer run 
showed that the ‘single 0 inversion program could not fit the 
histograms that were coming from the signal processor. The 
algorithm could not produce a fit because the lower signal bins 
were empty. The assumption made for all flights was that the 
threshold was set at 4 and therefore there would be counts in bin 
5 and above. For flight 13 the counts did not pick up until 
around bin 30. This left a hole in the particle cross-section 
distribution which the inversion algorithm could not fit with a 
power law, and spurious predictions resulted. 
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4.0 Conclusions 


The basis of the volume mode and single particle mode 
algorithms has been discussed. Programming and data handling 
techniques have been described. Some sources of systematic 
error in the previously processed flight data have been examined. 
Statistical uncertainty in the processed flight data has been 
evaluated. 

A rescaling of the single particle gain function as a 
result of a more exact handling of the calibration has resulted 
in a decrease by 2.9 orders of magnitude of the backscatter 
coefficient in the single particle mode. This leaves a differ- 
ential of about two orders of magnitude with the prediction of 
the volume mode method. Two possible explanations were offered 
for this difference: a) The single particle mode signal versus 
bin calibration for the processor is not correct, b) The particle 
density conditions during data collection did not conform to one 
particle per transit time. A further possibility is that the 
volume channel signal is either incorrect, or incorrectly inter- 
preted. 

The statistical errors for the two methods produced a 
maximum of ±.44 orders of magnitude uncertainty for the volume 
mode and +.54 orders of magnitude uncertainty for the single 
particle mode. The processed results do not agree within these 
errors . 

Studies of the effects of particle number fluctuations in 
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the processing of both modes show that, in order to reduce 
deviations of the local mean from a mean calculated with 10000 
particles, and to reduce fluctuations of the data from the local 
mean, care must be exercised regarding the number of particles 
processed. To keep these errors near .1 orders of magnitude, 
this number of particles must be 1000 for the single particle 
mode, and 3000 for the volume mode. This requires a very short 
flight length (<lkm) for particles with densities as low as one 
particle per 1000 cubic centimeters. 

In order to better understand the behavior of LDV systems 
operating in these two modes, it is suggested that laboratory 
measurements on aerosol standards be accomplished with the 
inflight processor. The calibration techniques which were 
interrupted by the flight schedule could be continued, or perhaps 
the input or exhaust from a Knollenberg counter could be used. A 
comparison of the two methods needs to be accomplished and their 
limits of applicability determined in the laboratory. The 
methods are independent assessments of aerosol distribution 
characteristics and, operating together, offer a unique means of 
converging on more precise measurements. 
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APPENDIX 


Program Descriptions 

Main Program - Sets up calibration data (except for gain versus 
area table); controls reading of data and number of 
particles/inversion; contains volume 6 computations and 
writes the output data file; computes flight length and 
single 3 after calling DATANAL. 

Subroutine GETDAT - Reads the flight data from disk as 192 words 
(4 bytes long) and converts to 384 words (unpacks the 
data). 

- Inputs - None 

- Outputs- BUF= 384 word array containing all data. 

DATA= 256 word array containing the counts in 
the signal processor for the single 
particle. 

Subroutine DATANAL - Sets up the gain versus area table and then 
calls OPTIM to do the actual inversion. 

- Inputs - NT= Total # of particle counts 

M = # of signal bin segments 

MS= # of cross-section segments 

- Outputs- WS= average cross-section 

WV= average focal volume cross-section area 

Subroutine OPTIM - does the inversion using three different fits 
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for the cross-section distribution: (1) power law (2) 

exponential (3) and log-normal. This routine iterates 
using each distribution type in turn to find the best 
match between the real single particle count histogram 
and the predicted count histogram. 

- Inputs - B= first power law exponent trial 

N= total # of particle counts 
M= # of signal bin segments 
MS= # of sigma (cross-section) segments 
NI= # of integration steps 

YY= array that contains scaled gain values for 
the "area" table 

ARR= array that contains areas for corresponding 
array YY 

- Outputs- WSS= predicted average cross-section 

WVS= predicted average focal volume cross- 
sectional area 

Subroutine C0EF2 - calculates the difference between two area 
values in an integration step. 

- Inputs - 1= index in SGL array (signal bin) 

J= index in SIG array (cross-section) 

K= integration step 
DJ= A sigma for index J 
M= # of cross-section segments 

- Outputs - C= area difference 


54 



Subroutine LPICK - generate a good starting point for the log- 
normal solution. 

- Inputs - SIG= cross-section array 

MS= # of cross-section segments 

- Outputs - BP= selected mean value of log 

DB= mean value of log increment selected 
ALP= selected standard deviation of log 
DALP= standard deviation of log increment selected 

Subroutine EPICK - generate a good starting point for the 
exponential solution. 

- Inputs - SIG= cross-section array 

BP1= selected power law fit from power law 
solution 

MS= # of cross-section segments 

- Outputs - BP= selected exponential power 

DP= exponential power increment 

Subroutine PROW - computes normalized probabilities for the power 
law, exponential and log-normal distribution depending 
upon the value of IP. 

- Inputs - BP= fit parameter for solution: 

• power law - exponent 

• exponential - exponential coefficient 

• log-normal - mean value of log 
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IP= l=power law 
’2=exponential 
3= log-normal 

SIP= selected cross-section for which probability 
will be calculated 

S0= (required only for log-normal) selected 
standard deviation of log value 
A3= (required only for log-normal) normalization 
factor for the selected log-normal 
parameters 

Function YLIN - linear interpolation function. 

- Inputs - N= # of data parts in array 

XX= dependent variable that results is derived 
to be interpolated for 
X= array of dependent variable values 
Y= array of independent variable values for 
which interpolations are performed 

- Outputs -YLIN= interpolated value 
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